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Abstract 

The possibility to improve the CERN to Gran Sasso neutrino beam performances 
for #13 searches is investigated. We show that by an appropriate optimization 
of the target and focusing optics of the present CNGS design, we can increase 
the flux of low energy neutrinos by about a factor 5 compared to the current 
r optimized focalisation. With the ICARUS 2.35 kton detector at LNGS and 
in case of negative result, this would allow to improve the limit to sin 2 2#i3 
by an order of magnitude better than the current limit of CHOOZ at Am 2 ~ 
3 x 1CT 3 eV 2 within 5 years of nominal CNGS running. This is by far the 
most sensitive setup of the currently approved long-baseline experiments and is 
competitive with the proposed JHF superbeam. 



1 Introduction 



The firmly established disappearance of muon neutrinos of cosmic ray origin [|l], |2| 
strongly points toward the existence of neutrino oscillations. 

The approved first generation long baseline (LBL) experiments — K2K 0, MI- 
NOS ICARUS || and OPERA || — will search for a conclusive and unambiguous 
signature of the oscillation mechanism. They will provide the first precise measure- 
ments of the parameters governing the main muon disappearance mechanism. In par- 
ticular, the CERN-NGS beam[]7], ^J, specifically optimized for tau appearance, aims to 
directly confirm the hints for neutrino flavor oscillation. 

The physics program of ICARUS will start with the installation of the 600 ton 
prototype at the Gran Sasso Laboratory (LNGS) and will allow the observation of 
atmospheric neutrinos, the detection of solar and supernovae neutrinos and the search 

1 On leave of absence from INFN Milano. 



for proton decay. An extension of the mass of argon is foreseen with the goal of reaching 
a total mass of about 3000 tons. 

In addition to the dominant —> v T oscillation, it is possible that a sub-leading 
transition involving electron-neutrinos occur as well. In the "standard interpretation" 
of the 3-neutrino mixing, the —> v e oscillations at the Am 2 pa 2.5 x 10 -3 eV 2 
indicated by atmospheric neutrinos is driven by the so-called 6*i 3 angle. Indeed, given 
the flavor eigenstates v^ia = e, /j,, r) related to the mass eigenstates u[{i = 1,2,3) 
where v a = U a iv[, the mixing matrix U is parameterized as: 

C12C13 S12C13 s^e ll 

U(0 12 , 13 , 9 23 , 8) = \ ~S 12 C 23 - C 12 S 13 S 23 e tS C 12 C 23 - S 12 Si3S 23 e lS C13S23 J (1) 

si 2 s 23 - ci 2 si 3 c 23 e lS -C12S23 - si2Si 3 C2 3 e J<5 Ci 3 c 23 

with Sij = sin and Cj 3 - = cos 9ij . 

The best sensitivity for this oscillation is expected for ICARUS at the CERN- 
NGS. Limited by the CNGS beam statistics at low energy, this search should allow to 
improve by roughly a factor 5 (see Ref.||) the CHOOZ|9[] limit on the #13 angle for 
Am 2 pa 3 x 10~ 3 eV 2 . Beyond this program, new methods will be required in order to 
improve significantly the sensitivity. 

At present, the only well established proposal in this direction is the JHF-Kamioka 
project In its first phase, 5 years of operation with the Super-K detector, it aims 



to a factor 20 improvement over the CHOOZ limit. 

In this paper, we investigate the possibility to perform 13 searches with a CERN 
to Gran Sasso beam optimized for Low Energy (L.E.) neutrino production. This is 
meant to be alternative to the CNGS r program, in order to be competitive with JHF. 

FLUKAjll], [12[] Monte Carlo simulations are employed to calculate neutrino yields 



with specific target and focusing configuration. 

Simulated neutrino fluxes are employed to derive exclusion plots in the Am 2 vs 
sin 2 (2^13) plane for the ICARUS detector. 



2 The neutrino energy range 

In order to maximize the probability of an oscillation, we must choose the energy of 
the neutrino E max and the baseline L such that 

L(km) 9, Tr9 . 71 

However, in order to observe the oscillation, we must at least see the maximum preceded 
by a minimum, given by 

E min (GeV) 
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Table 1: Neutrino energies E max and -E m in (see text for definition) corresponding to 
the maximum and minimum of the — > v x oscillation for the CERN-LNGS baseline 
and the Am 2 range indicated by Superkamiokande. 

At the CERN-LNGS baseline (730 Km), with Am 2 in the range indicated by Super 
Kamiokande (1 x 10~ 3 < Am 2 < 4 x 10~ 3 eV 2 ), the maxima and minima of the 
oscillation lie between 0.3 GeV and 2.4 GeV (See Table |l|). 

3 The present CERN-LNGS configuration 

The present CNGS design[[F| is optimized for v T appearance, thus for a relatively high- 
energy neutrino beam. The 400 GeV/c SPS beam will nominally deliver 4.5 x 10 19 
protons per year on a graphite target, made of spaced rods to reduce the re-interaction 
rate within the target. The two magnetic horns (horn and reflector) are tuned to focus 
35 and 50 GeV/c mesons, with an acceptance of the order of 30 mrad. The present 
shielding and collimator openings would not allow more than 100 mrad even in perfect 
focusing. The decay tunnel length is 1 km, and the baseline for neutrino oscillation is 
732 km. 

4 The CNGS L.E. option 

To improve particle yield at low energies, we re-designed the focusing system, and 
changed the target dimensions and shortened the effective decay tunnel length. The 
main differences wit the present (r) design are summarized in Table |2|. 

The L.E. target is still made of graphite, but the spaces have been eliminated. A 
more compact target ( 1 meter full length) allows for a better focusing and increases 
the low-energy yield as will be discussed in the next section. 

The decay tunnel length in the present calculations has been set at 350 meters. 
This corresponds to about two decays lengths for pions producing 2.5 GeV neutrinos. 

The neutrino energy of interest correspond to pions in the range 0.7-5.5 GeV. To 
focus these pions, we adopt a standard double-horn system (following the CNGS tra- 
dition, we call horn the first magnetic lens and reflector the second). Both have to 
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Target 

Material 

Total target length 
Number of rods 
Rod spacing 
Diameter of rods 

Horn 

Distance beginning of target-horn entrance 
Length 

Outer conductor radius 
Inner conductor max. radius 
Inner conductor min. radius 
Current 

Reflector 

Distance beginning of target-reflector entrance 
Length 

Outer conductor radius 
Inner conductor max. radius 
Inner conductor min. radius 
Current 

Decay tunnel 

Distance beginning of target-tunnel entrance 

Length 

Radius 



CNGS t 

Carbon 
2 m 
13 

first 8 with 9 cm dist. 
first 2 5 mm, then 4 mm 



320 cm 
6.65 m 
35.8 cm 
6.71 cm 
1.2 cm 
150kA 



43.4 m 
6.65 m 
55.8 cm 
28 cm 
7cm 
180kA 



100 m 
992 m 
122 cm 



CNGS L.E. 

Carbon 
1 m 
1 

none 
4mm 



25 cm 
4 m 

80 cm t 
11.06 cm 
0.2 cm 
300kA 



6.25 m 

4 m 

90 cm t 
23.6 cm 

5 cm 
150kA 



50 m 
350 m 
350 cm t 



Table 2: Parameter list for the present CNGS design and the "new" beam for low 
energy v . For the parameters flagged with a a full optimization has not been 
performed and possible improvements have not been studied yet. 

be placed near to or even around the target, to capture particles emitted at relatively 
large angles. The average transverse momentum of secondary particles is around 300 
MeV/c, with a sizable number of events up to 600 MeV/c, corresponding to 750 mrad 
for 1 GeV/c pions. These particles have to be bent before they travel too far away in 
radius, therefore the horn magnetic field has to be high enough. This also means that 
the horn should be shorter than the ones used to focus high energy beams, because 
the particles should not travel in the magnetic field for a distance longer than their 
curvature radius. 

We obtained good focusing capability with two four meter long horns. The horn 
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Table 3: Neutrino beam parameters for the CNGS baseline, with E v < 2.5 GeV. The 
t cases correspond to the present CNGS design for target, acceptance and focusing 
system. 



current has been set at 300kA, the reflector one at 150kA. The horn starts 25 cm 
after the target entrance face, the reflector starts just two meters after the horn end. 
Horn and reflector shapes has been computed to focus 2 GeV/c and 3 GeV/c particles 
respectively. We are aware that these (parabolic) horn shapes are derived in the ap- 
proximation of point-like source, that is not verified in the present case. However, the 
Monte Carlo calculations verified the good focusing capabilities of this system. The 
focusing efficiency in the range of interest is around 50%. 

5 Neutrino fluxes and rates 

With the standard CNGS parameters, the low-energy neutrino flux is low, as can be 
seen from the entries flagged by t in Table 0, even assuming perfect focusing: we expect 
0.9 CC events per kton per 10 19 pots for the tau focusing. Even if ideal focusing is 
assumed, the rate is only improved by a factor 2. 

The more compact target, and a wider acceptance, give about a factor five higher 
rate in the 0.-2.5 GeV range, or about 4.5 CC events per kton per 10 19 pots for the 
real focusing and up to 9.0 CC events per kton per 10 19 pots for the ideal focusing. 

While in both cases the focusing introduces an efficiency factor of about 50% with 
respect to the ideal focusing, our improvement comes from the ability to capture more 
wide angle, soft pions. The interplay of target configuration and angular acceptance 
can be seen clearly from Figure [1]: the low-energy part of the produced pion spectrum 
is enhanced due to the higher re-interaction probability in a compact target, at the 
expenses of the medium-high energy part, the one of interest for r appearance. This 
enhancement is small in the forward direction (100 mrad acceptance), but it grows at 
larger angles (up to 1 rad), where it reaches an average factor of 1.5 for pions in the 
1-5 GeV energy range. The real boost in low energy pion production comes however 
from the angular acceptance of the system, which accounts for a factor 3.3 in the case 
of a compact target and a factor 2.3 in the case of the CNGS target. The difference 
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in acceptances (1 rad in the low-energy option vs. 100 mrad for CNGS in perfect 
focusing) is therefore the dominant factor. 

Even after focusing, the L.E. option produces five times more low energy than 
the standard one. The difference can be appreciated by eye in Figure |2|, where we 
show the Up fluxes at Gran-Sasso for the two configurations. The electron neutrino 
contamination is of the same order of magnitude (around 1%) for the r and L.E. 
options. Electron neutrino spectra are plotted in Figure §. 

6 A near detector ? 

During the optimization for low energy pions, it came natural to reduce the "used" 
length of the decay tunnel from the available 1 km of the CNGS, in order to reduce 
the high energy neutrino tail. 

The reduction of the actual effective decay length could be accomplished by re- 
locating some of the graphite and iron blocks of the currently planned CNGS beam 
dump within the decay tunnel. We are aware that this operation is technically not as 
trivial as it might sound, however, we want to point out that it is not necessary to fill 
the entire decay tunnel but simply to stop the hadrons at the given point of the decay 
tunnel in order to achieve our wanted result. 

The reduced use of the decay tunnel leaves the room for an eventual "near" detector 
to be placed in the present beam dump position, that is at about 1 km from the target. 
This would allow us to monitor the beam in absence of oscillations and to predict the 
beam spectrum at the far position. 

As is well known however 0, ^, [T(|, the neutrino spectra at the near position are 
different from the "far" ones, mainly because they feel the finite size of the decay 
tunnel. We have simulated for the L.E. beam the fluxes as a function of v energy 
at 1 and 2 km from the target. The 1 km distance would correspond to the current 
"beam dump area" while the 2 km corresponds to the originally planned location of the 
TOSCA detector. We have rescaled our results according to the square of the baseline 
ratio, such that the spectra are directly comparable. As shown in Figures |] and |5|, the 
low energy part of the spectrum is enhanced with respect to the one at Gran Sasso, 
with differences up to 60% in the case of the very near detector. 

Therefore, the of v e beam contamination at Gran Sasso cannot be directly evaluated 
from the measurement at the near detector, but has to be "propagated" through Monte 
Carlo simulations. Hence, the need of the near detector is not fundamental, however, 
it would allow to cross-check the beam simulation in the region of no oscillations which 
can then be extrapolated at far distance. The resulting systematic error has not been 
estimated yet, however, we point out in next section that our results will be limited by 
statistics. 
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Table 4: Events from the CNGS L.E. beam, assuming Am|, = 3 x 10~ 3 eV 2 , 9 23 = 45°, 
5 years of operation and 2.35 kton fiducial mass. 
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Table 5: Events from the CNGS r beam, assuming Am 2 23 = 3 x 10~ 3 eV 2 , 6 23 = 45°, 5 
years of operation and 2.35 kton fiducial mass. 

7 Oscillated events and limits 

Figure |6| shows the exclusion plots for sin 2 2#i3 at the 90% C.L. for the CNGS r and 
our optimized CNGS L.E. options, compared with the CHOOZ limit, the SuperK 
allowed region and other proposed experiments. These contours have been derived 
from u e appearance, assuming 5 years of operation at the nominal CNGS intensity, 
and an ICARUS fiducial volume of 2.35 kton. For MINOS we assume an exposure of 



10 kton x year|14| and we quote the JHF limit according to the proposal of the OAB 
beam[]TU|. Three neutrino formalism and maximal v^-v T mixing have been assumed, i.e. 
6*23 = 45°, so that the usual "appearance factor" Piy^ — > u e ) = 0.5 x P(u e — > u^/u T ) 
has been taken into account. 

Only the intrinsic v e beam contamination has been considered here for the back- 
ground evaluation; other background sources, such as 7r° production in neutral current 
events, have been extensively studied in the past|| and found to be negligible given 
the excellent granularity of the ICARUS detector. They will nonetheless be precisely 
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re-addressed in a forthcoming paper ||15| . We however do not expect any significant 
changes in sensitivity introduced by other sources of backgrounds. 

The integrated number of background and oscillated v e events in the hypothesis 



Am 



3 x lO-W^ are listed in Table 



The oscillation maximum in the case 



of Am 2 = 3 x 10 _3 eV 2 , sin 2 20 13 = 0.1 is evident from the event spectra shown in 
Figure [8]. For comparison, we also quote the numbers for the CNGS r optimization in 
Table ||. 

The CNGS L.E. options represents an increase of roughly a factor 2 in sensitivity 
with respect to the CNGS r beam, within the SuperK allowed region. For Am 2 = 
2.5 x KT 3 eV 2 , we find: 



(sin 2 26 13 )cngs,t < 0.04 or 13 < 6° 
(sin 2 26 13 )cngs,l.e. < 0.02 or 9 13 < 4° 



(4) 
(5) 



The improvement over the CHOOZ limit is almost tenfoldH], and the performances 



are three times better than those foreseen by MINOS pl|: 

(sin 2 26 13 )chooz < 0.14 or 13 < 11° 
(sin 2 26 13 ) MINO s < 0-06 or 13 < 7° 



(6) 
(7) 



Finally, the JHF proposal with the OAB beam gives (TO 

(sin 2 20 13 ) jhf,oab < 0.006 or 9 13 < 2.2° (8) 

The comparison with the first phase (5 years) of the superbeam JHF is still slightly 
unfavourable, however we point out the possible time schedule and the probable SPS 
proton beam intensity upgrades. 



• The JHF-Kamioka experiment will start one or two years later than the CNGS 
beamQ; this allows to set at 7 years the CNGS L.E. data taking. 



As for proton intensity, it is expected that PS and SPS upgrades ||16|| will bring the 
accelerated intensity from 4.8 x 10 13 to 7 x 10 13 protons per cycle. This represents 
an increase of ~ 50% in the p.o.t. per year. This intensity upgrade is already 
taken into account in all the design specifications for the CNGS facility |T^]. We 
call this CNGS1.5. 



Accordingly, the sensitivity contour for 7 years data taking at CNGS1.5 nominal 
intensity is shown in Figure [7[ This gives for Am 2 = 2.5 x 10~ 3 eV 2 

(sin 2 29 13 )cngsi.5,l.e. < 0.015 or 13 < 3.5° (9) 



assume the starting dates of 2007 for CNGS and 2009 for JHF. 
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Since we are limited by the intrinsic v e background, the sensitivity scales with the 
square of the exposure and one would need to increase the intensity of the CNGS 
and/or the total mass of ICARUS by a factor 



in order to reach exactly the level of sensivity of the JHF superbeam. 

8 Conclusions 

The CERN-NGS has been originally optimized and coupled to high quality detectors 
in order to unambiguously give evidence for the — > u T flavor oscillation mechanism. 

However given the CERN financial situation, the CERN-NGS programme, approved 
in 1999, is now foreseen to start with full intensity in the year 2007. Given these 
delays, it is worth wondering if the priority of the program should be tau appearance 
("a confirmation of the oscillation mechanism") or electron appearance ("a potential 
discovery of a new flavor transition"). 

In this paper, we have studied a different optimization of the CNGS beam optics 
that is optimized for — > v e oscillation searches at the Am 2 indicated by the atmo- 
spheric neutrinos. We find that this beam coupled with the approved ICARUS T3000 
experiment at LNGS would offer great opportunities to find neutrino oscillations driven 
by 6*i3 on the same time scale as that of the proposed JHF superbeam. 

To reach a sensivity sin 2 2#i3 < 0.006, the total exposure should be increased by a 
factor 6 compared to what we have assumed here. This would require a big increase 
in beam intensity or a substantial increase of the liquid argon mass. The cost of 
multiplying the mass of ICARUS by a factor six is interestingly on the scale of the 
price of a new superbeam. 
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Figure 1: Pion production at 400 GeV, for the CNGS target within 100 mrad accep- 
tance, for a compact target ( the one of this study) within lOOmrad acceptance and 
within 1 rad acceptance 
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Figure 2: Muon neutrino fluxes at Gran Sasso, for the present CNGS design and for 
the new target/optics configuration. Whole spectra (top) and low-energy(bottom) are 
shown 
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Figure 3: Electron neutrino fluxes at Gran Sasso, for the present CNGS design and 
for the new target/optics configuration. Whole spectra (top) and low-energy(bottom) 
are shown 
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Figure 4: Muon neutrino fluxes at Gran Sasso and at "near" positions, one and two 
kilometers from the target. For direct comparisons, the flux for the near detectors have 
been rescaled according to the square of (target-near)/(target-Gran Sasso) distances. 
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Figure 5: Ratio of near/far muon neutrino fluxes. 
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MINOS, 10 kton-years 
ICARUS, 2,55 kton. 5 years CNGS, ^+4,6 
JHF+SK, 22.5 kton, 5 years 
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Figure 6: Expected sensitivity to sin 2 2#i3 mixing angle at nominal CNGS intensity, 
compared to existing results from CHOOZ§| and SuperK|l3| and expected results from 
MINOS |] and JHF-SK[|T§. 
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Figure 7: Comparison of expected sensitivity to sin 2 2#i3 mixing angle with an improved 
CNGSxl.5jl6|l and 7 years of running. 
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Figure 8: Comparison of v e CC spectra at Gran Sasso, in presence and in absence of 
oscillation. 
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